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Background: The lead region, Po, Pb, Hg, and Pt, shows up the presence of coexisting structures having different deformation 
and corresponding to different particle-hole configurations in the Shell Model language. 

Purpose: We intend to study the importance of configuration mixing in the understanding of the nuclear structure of even-even 
Po isotopes, where the shape coexistence phenomena are not clear enough. 

Method: We study in detail a long chain of polonium isotopes, 190 ~ 208 Po, using the interacting boson model with configuration 
mixing (IBM-CM). We fix the parameters of the Hamiltonians through a least-squares fit to the known energies and 
absolute B(E2) transition rates of states up to 3 MeV. 

Results: We obtained the IBM-CM Hamiltonians and we calculate excitation energies, B(E2)’s, electric quadrupole moments, 
nuclear radii and isotopic shifts, quadrupole shape invariants, wave functions, and deformations. 

Conclusions: We obtain a good agreement with the experimental data for all the studied observables and we conclude that 
shape coexistence phenomenon is hidden in Po isotopes, very much as in the case of the Pt isotopes. 

PACS numbers: 21.10.-k, 21.60.-n, 21.60.Fw 
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I. INTRODUCTION 

Shape coexistence has been observed by now all through the nuclear mass region, encompassing both light nuclei 
( 16 0 region) |lj], proceeding up to the region of very heavy nuclei in the Pb region Hi , and have been reviewed in a 
number of papers over a period spanning about 3 decades M- Recent advances in experimental methods to explore 
nuclei, removed far from the region of /3-stable nuclei, have opened up possibilities to explore the appearance and 
behavior of shape coexistence in series of isotopes and isotones Q- A wealth of new data^on both energy systematics, 
but, more important, on observables such as in-beam spectroscopy and lifetime data [§,ill, Coulomb excitation using 
inverse kinematics [10], direct reactions on unstable nuclei llj, radioactive decay modes at the limits of nuclear 
stability jl2j | , and break-up reactions jh'if have been at the basis of recognizing the rather universal appearance of 
shape coexisting phenomena Q. Moreover, measurement of the essential ground-state properties, such as masses & 
charge radii [151 ] and nuclear moments [161 ] . as well as the possibilities to study monopole E0 transitions in nuclei 
0 fl8j . helped completing the data basis in such a way as to confront theoretical modelling of nuclear structure 
properties in greater detail than before. 

From a theoretical side, present-day methods starting from both the nuclear shell model, or approaching the atomic 
nucleus using mean-field methods have resulted in developments of both new algorithms as well as making use of the 
increased computing possibilities (see 0-0 and references therein). The present status has evolved in a situation 
where the conditions for shape coexistence to occur are becoming understood. It looks like a balance between two 
opposing nuclear force components, i.e., on one side the stabilizing effect caused by the presence of closed shells (the 
monopole part), aiming at stabilizing the nucleus into a spherical shape, versus the low-multipole (mainly quadrupole) 
components redistributing protons and neutrons into a deformed shape is at the origin of the appearance of shape 
coexistence in a given mass region. Recent large-scale shell-model studies (LSSM), using diagonalization in a very large 
many-open shell basis in various mass regions [22[ or making use of an advanced Monte-Carlo shell-model approach 
(MCSM) [2lj| . have been carried out. Besides, the concept to start from deformed average potentials and calculating 
the total energy curves as a function of deformation has been explored, in particular, for nuclei in the Pb region 
[23M25I ] . However, recent studies flii . using a microscopic approach to determine the optimal mean fields, ev en g oing 
beyond when bringing in the nuclear dynamics, have given quantitative results - using both Skyrme forces [26l - |29[ and 
Gogny forces I30l - l34| . as well as making use of a relativistic mean-field approach I.'IT I f 01 - that are indicative of the 
above mechanism. Moreover, attempts have been made and are still improved to extract a Bohr Hamiltonian [4ll , |42| 
starting from a microscopic basis [43H45[ . 

From a microscopic shell-model point of view, the hope to treat on equal footing the large open neutron shell 
from N=126 down to and beyond the mid-shell N=104 region, jointly with the valence protons in the Pt, Hg, Po, 
and Rn nuclei, even including proton multi-particle multi-hole (mp-nh) excitations across the Z=82 shell closure, is 
beyond present computational possibilities. The truncation of the model space, however, by concentrating on nucleon 
pair modes (mainly 0 + and 2 + coupled pairs, to be treated as bosons within the interacting boson approximation 
(IBM) [0), has made calculations feasible, even including pair excitations across the Z=82 shell closure 0, 0] 
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FIG. 1: (Color online) Experimental energy level systematics for the Po isotopes. Only levels up to E x ~ 3.0 MeV are shown. 
The symbol 4 indicates the beginning of the region where extra levels of low-spin, indicated with spin, parity assignments 
between brackets, in the Nuclear Data Sheet references, start appearing. The symbol * for the 8 + level for mass A=202 
indicates that the energy is within an interval < 40 keV above the 6 + level. Dashed lines connect states that are supposed to 
have a similar structure. 


in the Pb region in a transparent way. More in particular, the Pb nuclei have been extensively studied giving rise 
to bands with varying collectivity depending on the nature of the excitations treated in the model space [49 54. 
More recently, detailed studies of the Pt nuclei have been carried out [55l - [6ll ] as well as for the Hg nuclei |62| - |64| 
in an attempt to describe the large amount of low-lying states and their E2 decay properties, explicitly including 
particle-hole excitations across the Z=82 shell closure. A novel mapping procedure to determine an algebraic IBM 
Hamiltonan has been proposed by Nomura et al. [65|, [66j. In contrast to the standard shell-model to boson model 
mapping method [67i. ]t38|. it has been shown that an IBM Hamiltonian can be determined, mapping a self-consistent 
mean-field total energy surface E(/3,7) (over the full /? — 7 plane) onto the corresponding IBM mean-field energy. Very 
recently, this method was extended to include intruder mp-nh configurations, with a detailed coverage of the Pt, Pb 
and Hg isotopes [69l472l | . 

Recently, a lot of new experimental results have become available in the Pb region, encompassing, besides the Pb 
isotopes, the nearby Hg, Pt and Po isotopes (see Section [XT] for references on the Po nuclei). In both the Pb and Hg 
isotopes, there is overwhelming information (excitation energies, B(E2) values, isotopic shifts, a-hindrance factors, 
etc.) by now that highlights the presence of shape coexistence, which is associated with the presence of prolate, oblate 
and spherical bands in the case of the Pb nuclei, and with the presence of a prolate deformed band together with a 
less deformed oblate structure, forming the yrast band in the case of Hg [S:] . Whereas the intruder bands are easily 
singled out for the Pb and Hg nuclei in which the excitation energies display the characteristic parabolic pattern 



































































































3 


with minimal excitation energy around the N=104 neutron mid-shell nucleus, this structure is not immediate in both 
the Pt and the Po nuclei. Therefore, in the present paper, we carry out an extensive study of the Po nuclei within 
the context of the interacting boson model, including 2p-2h excitations across the Z=82 proton closed shell, thereby 
extending the regular IBM model space, containing N bosons, with the intruder N+2 IBM model space, also taking 
into account the interaction between both subspaces, which is called the IBM-CM approach. 

The paper is organized as follows: in Section [Til we describe the experimental situation in the Po isotopes whereas 
in Section cm we present the various theoretical approaches that have been used in the literature to study the Po 
nuclei; in Section HVl we succinctly present the IBM-CM formalism as well as the fitting methodology used; here, we 
also discuss the main outcome of the calculations on energy spectra, electric quadrupole properties (B(E2) values, 
quadrupole moments), and its comparison with the available experimental data; in Section 0 we discuss the results 
on a-hindrance factors as well as on the isotopic shifts; Section [VI] is devoted to the description of nuclear deformation 
properties of Po nuclei as derived from the IBM-CM mean-field energy, from the study of the quadrupole shape 
invariants and from the study of the kinematic moments of inertia, characterizing the yrast band structure of Po 
nuclei. Moreover, we compare the present results with the nearby Pb, Hg and Pt isotopes. Finally, in Section fYIIl 
both the main conclusions as well as an outlook for further studies are presented. 


II. EXPERIMENTAL DATA: SITUATION IN THE PO NUCLEI 

The even-even Po nuclei span a large region of isotopes, starting with the lightest presently know 190 Po nucleus 
(N=106), passing through the end of the shell, 210 Po, all the way up to N=134 at 218 Po. 

Many experimental complementary methods have been used to disentangle the properties over such a large interval. 
These nuclei are extensively covered in the Nuclear Data Sheet reviews for A=190 [73|| . A=192 ,[74j, A=194 [75j, A=196 
A=198 a A=200 '78], A=202 0, A=204 [13, A=206 0, A=208 [0, and A=210 gg, and span the region 
we concentrate on in the present paper. 

The experimental information (up to the end of the 90’s), down to A=192 (N=108), was obtained mainly using 
early fusion-evaporation reactions, followed by in-beam y-ray spectroscopy, with, in particular, information on energy 
spectra systematics for A=192,194 [Mil, A=194 [13, A=196,198 0, A=196,198 0, A=198 0, and A=198,200 
[9l| . This particular part of the mass region has also been analyzed using a and /3 decay studies of mass-separated 
Rn an At nuclei, resulting in information about excited 0 + states in the isotopes with A=196-202 jH} and in en ergy 
spectra of the isotopes A=200,202 0], as well as making use of a-decay studies from the Rn nuclei for A=198 f94| . 
Alpha decay studies, moving from the Po nuclei into the Pb nuclei, have been a major fingerprint, in particular, in 
view of the shape coexisting structure in the daughter Pb nucle i. There exists an extensive set of results that has 
been accumulated during a period of about two decades [95lTl09j . 

It was observed that mass A=194 (see Fig. [T]) indicated a break in the energy scale observed for the heavier masses, 
where the energy of the level is typically of the order of ~ 600 — 650 keV, going down to an energy of about ~ 300 
keV, therefore, dropping by a factor two. More recenty, experimental studies spanning the 1999-2009 period have 
been performed at the JYFL at the University of Jyvsk yla using the recoil-decay-tagging (RDT) technique, studying 
prompt and delayed y-rays down to mass A=190 [110;, 111]. Moreover, lifetimes could be derived makin g use of the 
RDT technique in recoil distance Doppler-shift experiments for masses A=194 [HD, tHH) and A=196 |l!4] . The 
experimental situation for the Po isotones down to A=190 has been summarized by Julin et al. [|]. Very recently, 
Coulomb excitation experiments have been performed using post-accelerated 196 > 19 8 , 200 , 202 p o beams at REX-ISOLDE, 
resulting in an important set of specific reduce d E2 matrix elements, connecting the 0]/’ ground state with the excited 
2+2 states (as well as for other combinations) lll5j . Charge radii have been measured for both the odd-A A=191-211 
Po nuclei ei, m3 as well as for the even-even A=192-210, 216, and 218 Po nuclei [l!8j |. extending the ground-state 
information in a major way. 


III. THEORETICAL APPROACHES: THE EVOLVING SITUATION IN THE PO NUCLEI 

The experimental situation, as discussed before, hints for a particular change in the structure moving down in mas s 
number A. The energy scale, which starts at the N=126 closed neutron shell, set by the seniority scheme EMU, 
with a first excited 2 + state at 1181 keV, quickly sets in a rather constant value of the 2 + energy at ~ 650 keV. 
However, over the span of a few mass units one notices (i) a quick drop of the first excited 0 + from mass A=202 
(1758 keV), down to mass A=196 (558 keV), and (ii) starting at mass A=196, the 2+ energy (at 463 keV) drops 
considerably setting into a new energy scale for the lower masses (known only down to A=190 at present), of ~ 250 
keV. 
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In the early works of [87j, [89j a description using a particle-core coupling model approach was proposed, resulting 
from coupling the 2 proton particles with a vibrational core sy stem. W ith the experimental information on masses 
A=196, and later A=194, a systematic study was carried out in jl 22 l fl23 1] for the mass region A=210 down to A=194. 
Here, it was pointed out that the quick drop of, in particular, the first excited 0^ state, arose be caus e of the need 
of using an enhanced particle-core coupling strength. This point was later stressed by Oros et al. |l24| who showed 
that a consistent set of coupling strengths was unable to describe the observed properties, in particular, below mass 
A=200. 

Calculations using a deformed average field |l25| indicated that in studying the energy surfaces, while at mass 
A=196, a soft spherical result was obtained, at mass A=192, an oblate minimum appeared, becoming the lowest 
minimum at a value of e =- 0.2 and an almost degenerate situation at A=190 with both oblate and prolate minima 
with both e values of ~ 0.2. Energy surface calculations, covering the /3 — 7 plane were carried out using a def ormed 
Woods-Saxon potential, resulting in the presence of three minima in the mass region A=190 to A=186 1241 ]. The 
above results were confirmed later by Smirnova et al. |l26| . this time using a self-consistent HFB approach using 
the SLy4 Skyrme force [2f|, indicating rather wide minima in A=196-194, a lowest oblate minimum at A=192 and 
A=190, turning into prolate at A=188. Trigger ed b y the new data on masses A=196 and A=194, more detailed 
beyond mean-field studies were carried out [l 12Ml 14 1 highlighting a detailed comparison on energy spectra and E2 
properties, indicating the interplay of an oblate and a more spherical structure (vibrational). In their analysis, the 
authors point towards a rather pure intruder character of the whole yrast band, including the ground state. In our 
analysis (see Section HVEl) , we come to the conclusion that the ground state exhibits a rather m ixed character instead 
of being of pure intruder nature. More recently, a full study of even-even nuclei in the Pb region ll 27i . showing specific 
results for the Po nuclei is presented. It is interesting to point out that high-spin isomers have been studied in the Po 
too, m aking use of a deformed Woods-Saxon potential, showing the effects of deformation on the specific excitation 


energy 


1128, 12S 


Most of the mean-field studies point towards the existence of rather complex energy surfaces with the presence of 
several minima, although in many cases separated by small barriers. Anyhow, a common denominator is the presence 
of a regular configuration, slightly deformed (oblate or gamma unstable) or spherical, coexisting with an intruder 
configuration of prolate nature, corresponding with a larger deformation as compared with the regular configuration. 
The lack of full calculations in the /? —7 plane makes difficult to formulate more precise conclusion about the particular 
shape of the coexisting minima. 

Within the framework of the shell-model and allowing for both the full neutron open shell, covering N=126 to 
N=82, as well as allowing proton mp-mh excitations across the Z=82 closed shell, the calculations are unfeasible 
nowadays. Therefore a truncated approach can be used starting from the standard IBM allowing for the presence 
of extra pairs. This method, called IBM-CM was proposed by Duval and Barrett [13, EH and has been used in the 
study of shape coexistence in various mass region. Some early studies were carried out within the idea of a possible 
symmetry to be used within an extended version of the IBM, includi ng p article-hole pairs (namely intruder I-spin 
Il30lll3l| b with specific applications to the Po nuclei (references jl24l . il32l ]l in which the coupling between U(5) and 
SU(3) symmetries were explored and compared to the then existing data (mass A=200 down to A=192)). Besides, a 
different symmetry, i.e., F-spin symmetry [[46], was proposed to relate energy spectra in nuclei with different number 
of protons, N w , and neutrons pairs, N„, outside of the nearest closed shells, ho wever keeping the sum F=(N W + 
N„)/2 constant. An application to the Pb region was carried out by Barrett et al. Il33| . The data obtained recently 


|134Hl40f from detailed studies of the Hg nuclei on energy spectra and electric quadrupole properties (B(E2) values, 
Q-moments), however, do indicate that the more simple idea of I-spin is not holding so well. 


IV. THE INTERACTING BOSON MODEL WITH CONFIGURATION MIXING FORMALISM 

A. The formalism 

The IBM-CM is an extension of the original IBM allowing to treat simultaneously several boson configurations 
which correspond to different particle-hole (p-h) shell-model excitations EH- In our case, the model space includes 
the regular proton 2p configurations and a number of valence neutrons outside of the Z=82, N=126 closed shells 
(corresponding to the standard IBM treatment for the Po even-even nuclei) as well as the proton 2h-4p configurations 
and the same number of valence neutrons corresponding to a [N] ® [N + 2] boson space (N being the number of active 
protons, counting both proton holes and particles, plus the number of valence neutrons outside the Z=82,N=126 
closed shells, divided by 2 as the boson number). Consequently, the Hamiltonian for two configuration mixing can be 
written as 


H = PIh^Pm 


+ P, 


iV+2 


(C q + f 2 + A n+2 )Pn+2 +V z 


iV,iV+2 


(i) 
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where Pn and Pjv +2 are projection operators onto the [TV] and the [TV+2] boson spaces, respectively, P^ Ar+2 describes 
the mixing between the [TV] and the [TV + 2] boson subspaces, and 

#ecqf = £ i fl d + k[L ■ L + KiQ(Xi) ' Q(Xi), ( 2 ) 

is a restricted IBM Hamiltonian called extended consistent-Q Hamiltonian (ECQF) jl4ll . Ill2j | with i = TV, TV + 2, fid 
the d boson number operator, 

L^iPxdp, ( 3 ) 

the angular momentum operator, and 

QuiXi) = [s t xT + d t x s]^ 2) + Xi[d ] x d}^> , (4) 

the quadrupole operator. This approach has been proven to be a good approximation in several recent papers on Pt 
[57l m an d Hg isotopes [sU. 

The parameter A N+2 can be associated with the energy needed to excite two proton particles across the Z=82 shel l 
gap, giving rise to 2 p- 2 h excitations, corrected for the pairing interaction gain and including monopole effects EMU. 
The operator describes the mixing between the TV and the TV + 2 configurations and is defined as 

yN ^+2 _ W ^’N+ 2 ( g t X gt + s X s) + W^ r,Ar+ 2 ((i^ xP+dxdr. ( 5 ) 

The E2 transition operator for two-configuration mixing is subsequently defined as 

T(E2) I1 = J2 ei p iQ»(Xi)Pi, ( 6 ) 

i=N,N+2 


where the (i = N, N + 2) are the effective boson charges and Q^(xi) the quadrupole operator defined in equation 

©. 

In section IIV Bl we present the methods used in order to determine the parameters appearing in the IBM-CM 
Hamiltonian as well as in the T{E2) operator. 

The wave function, within the IBM-CM, can be described as 

*(k,JM) = ^a a fc (J;TV)V»((sd)f;JM) 

i 

+ ^ fe (J;IV + 2)^((sd)f + 2 ;JM) , (7) 

3 

where k, i, and j are rank numbers. The weight of the wave function contained within the [TV]-boson subspace, can 
then be defined as the sum of the squared amplitudes w k {J , TV) = JT | a k ( J; TV) | 2 . Likewise, one obtains the content 
in the [TV + 2]-boson subspace. 


B. The fitting procedure: energy spectra and absolute B(E2) reduced transition probabilities 

Here, we present the way in which the parameters of the Hamiltonian (JT|), ([2]). (U), and J5|) and the effective charges 
in the T(E 2) transition operator (© have been determined. We study the range 190 Po to 208 Po, thereby, covering 
almost the whole second half of the neutron shell N=82-126. 

TABLE I: Energy levels, characterized by ,/,f, included in the energy fit, if known, and the assigned a values in keV. 


Error (keV) 

States 

cr = 0.1 

2 1 

(7—1 

4+,0+,2+ 

a = 10 

2^", 3f, 4^", 6j 1 ", 8^" 

a = 100 

6 2 + 
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In the fitting procedure carried out here, we try to obtain the best possible agreement with the experimental data 
including both the excitation energies and the B(E2) reduced transition probabilities. Using the expression of the 
IBM-CM Hamiltonian, as given in equation 0, and of the E2 operator, as given in equation ([6|), in the most general 
case 13 parameters show up. We impose as a constraint to obtain parameters that change smoothly in passing from 
isotope to isotope. Note also that we constrained \n = 0 and n' N+2 = 0. We have explored in detail the validity 
of this assumption and we have found very little improvement in the value of y 2 (see Eq. (|8]>) when releasing those 
parameters. On the other hand, we have kept the value that describes the energy needed to create an extra particle- 
hole pair (2 extra bosons) constant, i.e., A N+2 = 2400 keV, and have also put the constraint of keeping the mixing 
strengths constant too, i.e., w^’ N+2 = w 2 ' N+2 = 30 keV for all the Po isotopes. We also have to determine for each 
isotope the effective charges of the E 2 operator. This finally leads to 8 parameters to be varied in each nucleus. 

To determine the value of A w+2 we considered as a reference the already known values for Pt and Hg, which are 
2800 keV and 3480 keV, respectively. On one hand, we expect to get a value closer to Pt than to Hg because the 
intruder states for Po are supposed to be very close to the regular ones and even become the ground state around 
the mid-shell. That constrains our value to be not much higher than 2800 keV. On the other hand, we know that the 
heavier isotopes, near the closed shell, should have an energy gap between the regular and the intruder states which 
is equal, at maximum, to A w+2 (see Fig. [2] on correlation energy) and, moreover, experimentally the intruder states 
should appear above 2000 — 3000 keV. 

To determine the value of the mixing strengths, we considered that the corresponding value for the Pt nuclei was 
fixed to 50 keV [57|, while for the Pb, to a smaller strength of 18 keV [52, 53j, and for Hg to 20 keV. We performed 
a set of exploratory calculations between 20 — 30 keV and found that the best overall agreement corresponds to 
w^’ N+2 = w 2 ’ N+ 2 = 30 keV, although the agreement is quite similar over the whole range. 


TABLE II: Hamiltonian and T(E2) parameters resulting from the present study. All quantities have the dimension of energy 
(given in units of keV), except xn +2 which is dimensionless and ejv and ejv +2 which are given in units \AV.u. The remaining 
parameters of the Hamiltonian, i.e., xn, and k' n+2 are equal to zero, except for A^"*" 2 = 2400 keV and w^’ N+2 = w 2 ’ +2 = 30 
keV. 


Nucleus 

£N 

Kn 

k' n 

EN+2 

KN+2 

XN+2 

ei v 

ejv +2 

190 Po 

712.6 

-29.41 

6.34 

285.43 

-28.56 

0.22 

2.88“ 

1.86“ 

192 Po 

731.8 

-24.20 

0.43 

402.67 

-30.09 

0.21 

2.88“ 

1.86“ 

194 Po 

800.3 

-15.12 

-16.00 

518.23 

-28.92 

0.41 

2.88 

1.86 

196 Po 

845.0 

-22.76 

-10.96 

373.16 

-31.02 

0.15 

1.86 

1.86 

198 Po 

982.3 

-28.00 

-25.74 

854.82 

-34.51 

1.33 

2.05 

1.10 

200p o 

955.8 

-26.05 

-25.00 

843.48 

-34.55 

1.31 

2.01 

1.10 6 

202 Po 

942.1 

-29.06 

-28.46 

236.30 

-20.63 

1.09 

2.28 

1.10 6 

204 Po 

810.4 

-0.00 

-21.09 

100.00 

-5.00 

0.50 

2.28 c 

1.10 6 

206 Po 

717.1 

-0.00 

-4.30 

100.00 

-5.00 

0.50 

2.28 c 

1.10 6 

208 Po 

643.4 

-0.00 

6.52 

100.00 

-5.00 

0.50 

2.28 c 

1.10 6 


“The effective charges have been taken the same as the corresponding values obtained for 194 Po. 

“ejv +2 corresponding to 198 Po. 
c ej v corresponding to 202 Po. 

The x' 2 test is used in the fitting procedure in order to extract the optimal solution. The y 2 function is defined in 
the standard way as 

2 1 ^ (Xj(data) — Xj(IBM)) 2 

X Ndata - N p ar of 

where Ndata is the number of experimental data, N par is the number of parameters used in the IBM fit, Xi(data ) 
describes the experimental excitation energy of a given excited state (or an experimental B(E2) value), Xi(IBM) 
denotes the corresponding calculated IBM-CM value, and ct,; is an error (theoretical) a ssign ed to each Xi(data) point. 
We minimize the y 2 function for each isotope separately using the package MINUIT |l44l | which allows to minimize 
any multi-variable function. 

In some of the lighter Po isotopes, due to the small number of experimental data, the values of some of the free 
parameters could not be fixed unambiguously using the above fitting procedure. Moreover, for the heavier isotopes 
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(A>202), that part of the Hamiltonian corresponding to the intruder states is fixed such as to guarantee that those 
states appear well above the regular ones, that is, above 2 MeV. In some cases due to the lack of experimental data 
the effective charges could not be determined. 

As input values, we have used the excitation energies of the levels presented in Table Q] In this table we also give 
the corresponding a values. We stress that the a values do not correspond to experimental error bars, but they are 
related with the expected accuracy of the IBM-CM calculation to reproduce a particular experimental data point. 
Thus, they act as a guide so that a given calculated level converges towards the corresponding experimental level. The 
a (0.1 keV) value for the 2+ state guarantees the exact reproduction of this experimental most important excitation 
energy, i.e., the whole energy spectrum is normalized to this experimental energy. The states 4+,0 \ and 2^ are 
considered as the most important ones to be reproduced (a = 1 keV). The group of states 2(j", 3+, 4^", 6+, and 8+ 
(a = 10 keV) and 6^ (a = 100 keV) should also be well reproduced by the calculation to guarantee a correct moment 
of inertia for the yrast band and the structure of the 0^ band. Note that we only considered states in the fit with 
angular momentum and parity unambiguously determined. 

In the case of the E2 transitions rates, we have used the available experimental data involving the states presented 
in Table HI restricted to those E2 transitions for which absolute B(E2) values are known, except if serious hints that 
the states involved present non-collective degrees of freedom exist. Additionally, we have taken a value of a that 
corresponds to 10% of the B(E2) values or to the experimental error bar if larger, except for the transition 2+ —»• 0+ 
where a smaller value of a (0.1 W.u.) was taken, thereby normalizing, in most of cases, our calculated values to the 
experimental B(E2;2+ —> Of) value. 

This has resulted in the values of the parameters for the IBM-CM Hamiltonian, as given in Table [TT] In the case 
of 19 °- 192 p 0 anc j 200—208p 0 ^ ^ ie va j ue G f the effective charges, or part of them, cannot be determined because not a 
single absolute B(E2) value is known or y 2 is insensitive to their values. However, for completeness we have taken 
the effective charges of 194 Po for 19 °- 192 p 0; ejv +2 of 198 Po for 200 - 208 p Oj and e n of 202 Po for 204 - 208 p o . 



FIG. 2: (Color online) Absolute energy of the lowest unperturbed regular and intruder 0) 1 " states for 190 208 Po. The arrows 
correspond to the correlation energies in the N and N+2 subspaces (see also the text for a more detailed discussion). 

















C. Correlation energy in the configuration mixing approach 


Intruder states are expected to appear, in principle, at an excitation energy well above the corresponding regular 
ones with identical angular momentum. The reason is that these configurations are related to the creation of a 2p-2h 
excitation across the Z=82 closed shell. In the case of Po this energy would correspond to A iv+2 = 2400 keV. However, 
according to reference [ 13 , the single-particle energy cost has to be corrected because of the strong pairing energy 
gain when forming two extra 0 + coupled (particle and hole) pairs, by the quadrupole energy gain when opening up 
the proton shell, as well as by the monopole correction caused by a change in the single-particle energy gap at Z=82 
as a function of the neutron number. In particular, around the midshell point at N=104, where the number of active 
nucleons becomes maximal, the energy gain due to the strong correlation energy is such that the energy of the intruder 
configurations becomes close to the energy of the regular ones. It may even be that a crossing results making the 
intruder state forming a ground state. For instance, in the case of Pt isotopes, the nuclei around the midshell point 
at N=104 exhibit a ground state of intruder nature [57i, [58f. 
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FIG. 3: (Color online) Experimental excitation energies (up to E x ss 3.0 MeV) (panel (a)) and the theoretical results (panel 
(b)), obtained from the IBM-CM. 

The experimental states can present a strong mixing because of the interaction between both families of configura¬ 
tions (regular and intruder). Therefore, it is not simple in many cases to find out which configurations are dominant 
in the ground-state wavefunction. In order to take advantage of the IBM-CM calculations, we calculate explicitly the 
“absolute” energy of the lowest 0 + state belonging to both the regular [N] and [N+2] intruder configuration spaces, 
turning off the interaction among the two families. We choose as the reference energy, the energy of the regular 
Hamiltonian describing the [N] space for the 210 Po nucleus. This nucleus only consists of a single boson, with the 
s-boson as the lowest energy state, resulting in the zero-energy line (horizontal dashed red line in Fig. [5]). The energy 
of the lowest 0 + state in the regular configuration space [N], E(0+,N), is lowered with respect to the reference energy 
because of the correlation energy and is described by the wave function 4>(0) l ")^ s (see also expression ©). The 
lowering depends on the number of bosons. On the other hand, the energy of the lowest 0 + state in the intruder 
configuration space [N+2], E(0f~ ,N+2), is described by the wave function 4>(0^))v +2 ( see a ls° expression (fTOl) I and 
appears at the energy corresponding to A^" 1 " 2 . This energy will subsequently be lowered by its specific correlation 
energy too. In most cases, the regular configuration with N bosons corresponds to a spherical or slightly deformed 
shape, while the intruder ones, with N+2 bosons, to a more deformed shape. Therefore, the energy gain for the lowest 
intruder state, described by the wave function 'F(C)i~))v+ 2 > uses t° be larger than for the lowest regular state, described 
by the wave function T(0f ))^ 9 . The relative position of these lowest regular and intruder states is plotted in Fig. [2] 
Here, it can be clearly appreciated how the energies of both configurations can become very close, depending on the 
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balance between the off-set, A Ar+2 , and the difference in the correlation energy E(N + 2) corr — E(N) corr . 

One notices how the energies become really close near midshell (N=106, A=190), where the number of active 
nucleons is maximal, showing that the lowest intruder state can determine the character of the lowest 0+ state. From 
A=196, and moving towards the heavier masses, the energy difference starts to increase and both states appear to 
be very well separated. Note that once the isotopes approach the end of the shell, at N=126, the number of active 
bosons is drastically reduced and, therefore, the correlation energy is reduced dramatically. In particular, the energy 
of the regular state reaches the reference energy, which means that it corresponds to essentially a spherical shape, 
while the energy of the intruder state approaches the energy A w+2 quite closely. Consequently, at the end (or at the 
beginning) of the shell, the maximum energy difference between both states will correspond to A N+2 . The existence 
of a ceiling for the energy difference modulates the parabolic behavior of the energy systematics of the intruder states, 
transforming it into a flat shape and, therefore, leading to an energy for the intruder states lower than experimentally 
observed. Clearly, this is a deficiency of the IBM-CM calculations near the shell closure. 



FIG. 4: (Color online) Comparison of the absolute B(E2) reduced transition probabilities along the yrast band, given in W.u. 
Panel (a) corresponds to known experimental data and panel (b) to the theoretical IBM-CM results. 


D. Detailed comparison between the experimental data and the IBM-CM results: energy spectra and 

electric quadrupole properties 

In this subsection, we compare in detail the experimental energy spectra with the theoretical ones up to an excitation 
energy below E x ps 3.0 MeV. In Fig. [3jt we plot the experimental data set, while in Fig. [3jo the theoretical values are 
depicted. Note that the energy corresponding to the 2f state perfectly agrees with the theoretical one because this 
level was used to normalize the energies, in other words, we used in the fitting procedure for this level a very small 
value for a (0.1 keV) in order to exactly reproduce the 2^ experimental energy. Comparing the energy spectra, one 
can distinguish two regions, A < 200 and A > 200, with an overall better agreement in the first region as compared 
to the second. A most probable reason for this is the fact that in the lightest isotopes the number of active bosons 
is substantially larger and the collectivity is enhanced. As we discuss later, in the heavier Po isotopes, non-collective 
broken-pair states (with a v=2 seniority character) do appear in the low-energy part of the energy spectrum. Those 
states are not well reproduced within the IBM space because of the restriction to a set of interacting s and d bosons. 

One observes in Fig.[3]how the compression of the energy spectrum, with decreasing mass number, is well reproduced. 
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A A 

FIG. 5: (Color online) Comparison of the few non-yrast intraband absolute B(E2) reduced transition probabilities, given in 
W.u. Panel (a) corresponds to the few known experimental data, panel (b) to the theoretical IBM-CM results. 


The same happens for the systematic dropping of the energy of the 0^ states, and the energy of the yrast band, more 
generally. Note that the agreement is better for those states with low and even angular momenta. 

Next, we carry out a comparison for the B(E2) values, which is a much more stringent test than the excitation energy, 
because these numbers are highly dependent on the detailed structure of the wave function. Although the existing 
experimental information on B(E2) values is rather scarce, recent new re sults from Coulomb excitation experiments at 
REX-ISOLDE, in particular for the 196 - 198 p 0 isotopes, have appeared ina, improving the experimental knowledge 
of this mass region. 

In Figs. |4] and [5] we compare the B(E2) reduced transition probabilities, while in Fig. [6] we compare the electric 
quadrupole moments. We also present a more detailed comparison on B(E2) values in Table Hill 

In Fig. [3] we present the intraband B(E2) values along the yrast band, for which the most complete experimental 
information exists. In a similar way as for the excitation energies, we used the B(E 2; 2+ —> 0)*") reduced transition 
probability to normalize the theoretical results. Note that for some nuclei, where experimental information is scarce, 
we have used the corresponding values of the neighboring isotopes to determine the effective charges (as explained 
in Sec. IIV I ill . On the experimental side, the value of B(E 2; 2^ —>• Of) is systematically dropping, which denotes a 
reduction of the collectivity of the states when approaching the end of the neutron shell. However, this is not the 
case for B(E2\A J [ —> 2+) in 198 Po, where an unexpected large value is observed. On the other hand, the theoretical 
results show a continous drop of all the intraband transitions, resulting from the reduction of the collectivity of the 
states (in line with a decreasing number of active bosons). 

In Fig. [5] we study the systematics of some interband transitions involving the states and 2)*~ 2 . Because of 
the lack of experimental data and the large error bars, it is difficult to extract any trend. Regarding the theoretical 
results, once more, one notices an overall reduction of the B(E2) values. We point out that for A > 196, the observed 
B(E2) values are of the same order as the observed B(E2) values in the yrast band. This represents a hint about the 
changing mixing character of the states as a function of mass number A (see Fig. [12] as illustration of the changing 
structure of the wave function). 

A different way to extract information on the changing character of the wave functions for the lowest two 2 + 
states is comparing the spectroscopic quadrupole moments, shown in Fig. [6] Here, we compare the theoretical and 
the experimental values of the quadrupole moment for the states 2^ and 2^~. Moreover, we present the values 
corresponding to the unperturbed lowest-lying three 2 + states. One notices that the 2+ state corresponds to a well 
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TABLE III: Comparison of the experimental absolute B(E2) values (given in units of W.u.) with the IBM-CM Hamiltonian 
results. Data are taken from the Nuclear Data Sheets [734781 ]. complemented with references presented in section El 


Isotope 

Transition 

Experiment 

IBM-CM 

194 Po 

2+^0+ 

90(20) 

90 


4+^2+ 

120(40) 

156 

196 Po 

2+^0+ 

47(6) 

47 


4+^2+ 

103(15) 

92 


6+^4+ 

130(60) 

125 


2+ - 
Z 2 

+ 0+ 

5.7(10) a 

2.1 


2+ _ 
Z 2 

->2+ 

130(25) a 

36 

198 Po 

2+^0+ 

39(9) 

39 


4+^2+ 

177(53) 

61 


8+^6+ 

2.0(1) 6 ’ c 

72 


oj- 
2 1 - 

+ 2+ 

->o+ 

285(11®°)“ 

vs (Uir 

7 

1.7 


2+ ->• 2+ 

229(147)“ 

13 


2+ - 
Z 2 

->o+ 

42(56)“ 

33 

200 po 

2+^0+ 

31(2) 

31 


8+^6+ 

9.4(5) 6 

52 

202 po 

2+^0+ 

32(8) 

32 


“Data taken from Ref. Ill5l . 

“Experimental data not included in the fit. 
“Data taken from Ref. l9lll . 


deformed and oblate shape for the lightest isotopes, which is smoothly changing into a rather spherical shape for 
A=202 and onwards. On the other hand, the 2J state corresponds to a prolate shape for the lightest isotopes, 
changing into an oblate shape for A=196-200, and into a spherical shape for A=202 and the still heavier isotopes. 
When comparing with the unperturbed values, one can see how the 2+ state is built up as a mixture of the first two 
unperturbed 2+ states. However, in order to understand the second 2^ state, we have to resort to a mixture of the 
second and third unperturbed 2 + states mainly. This figure presents a clear proof of the changing character for the 
first two 2 + states as a function of mass number. This issue will be discussed in section ttV El in a more quantitative 
way. 

In Figs. [7] and [5] we present the experimental and theoretical energy spectra (up to E x ss 2.5 MeV) for masses A=190- 
198, which is the region where the coexistence should be more evident. We include in the comparison the known 
absolute B(E2) values. One observes a rather distinct rotational structure (along the yrast band) for A=190-194, 
followed by a structure that is changing into a more vibrational behavior for masses A=196-198. 

The Po energy spectra show some other most interesting structural changes starting at the neutron N=126 shell 
closure, moving towards the neutron-deficient region (decreasing towards mid-shell at A=188) (see Figs. [Q and [3]). 
One starts with the typical 7r(l/ig/ 2 ) 2 two-particle spectrum governed by the seniority characteristics of the two-body 
proton-proton interaction. An energy gap of ~ 0.8 MeV is observed from the 8+ to the 8^" state energy difference (see 
Fig. hd, which most probably results from a Tr(lhg/ 2 2f 7 / 2 ) configuration, as a measure of the single-particle energy gap. 
Next, one observes the remaining characteristics of a senority v=2 spectrum, albeit with an immediate drop of the 2 7 
state, when removing two neutrons from the closed N=126 shell, which then keeps a remarkably constant energy down 
to mass A=200 (N=116). The energy of the 4+ slighty drops but also remains at an approximate constant energy. 
The higher spin members 6+, 8+ are shifted upwards in a rather regular way down to mass A=198 (N=114), before 
first the 6)*" starts dropping systematically down in excitation energy at A=196 (N=112), followed by the 8+ dropping 
too from A=194 (N=110) (see Figs. [Hand 0 . Therefore, it looks like the region 208 > A > 198 is reminiscent of a 
vibrational pattern where the specific two-nucleon properties of the high-spin proton pair remain rather well intact 
down to the lower value at A=198. An interesting test to gain deeper insight in this region, and the connection to 
the yet lower mass region 196 > A > 190 where a dominant rotational-like and more collective pattern is showing 
up, can be derived from a study of, e.g., the B(E2;8) 1 " —> 6+) value, as well as, o f the g-factor and the quadrupole 
moment of the high-spin 8+ state (see Table Hill and Fig.lOj. Experiments |l45Lll4q have shown that the g-factor stays 














12 



FIG. 6: (Color online) IBM-CM values of the quadrupole moments for the 2+ and 2\ states in 190 ~ 208 Po. The quadrupole 
moments are given in units e b. The dash-dotted and dotted lines indicate the quadrupole moments when the mixing Hamiltonian 
is switched off, corresponding to the unperturbed 2+, 2^, and 2^" states. 


remarkably constant going down to A=198, at a value consistent with the 7r(l/ig/2) 2 8 + configuration (for a pure v=2 
configuration, as a function of n, the number of holes in the 7r(l/ig/ 2 ) orbital, the g-factor remains constant), indicating 
a minor influence of collective admixtures into this state. Maj et al. [9l| have analyzed the B(E2;8] I ~ —I 6+) value as 
well as the 8+ quadrupole moment in order to extract an effective charge (see e.g., Fig. 9 in [91(). The influence of 
particle-core coupling on the quadrupole moment of the 8+ state has been stu died for the Po nuclei with A=202 and 
204 indicating an increase in the quadrupole moment as low as A =200 |l47t |. The effective charge shows a steady 
increase of ~ 2 from A=210 down to mass A=200, from which a sharp drop is observed moving towards A=198. This 
seems to support the idea that, in particular for the high-spin members with a seniority v=2 character, those states 
only couple weakly with collective quadrupole excitations. Moreover, these results also corroborate the fact that from 
A=196 and onwards with decreasing A value, large collective compon ents shou ld appear in the 6^" state, even though 
the g-factor of the 8+ remains constant between A=200 and A=198 [l45l . Il4tj| . In view of the discrepancies between 
the IBM-CM calculations and the experimental data concerning the g-factor and the B(E2;8] I ~ —> 6) 1 ") value, those 
states should have a very weak coupling with the quadrupole collective excitations, and genuinly correspond to almost 
pure v=2 seniority configurations (see Fig. El). 

Finally, it is worth to note a recent beyond-mean-field calculation performed for the lead region by Yao et al. [ 1271] 
and, in particular, with results for the even-even Po isotopes (A=186-204). They have calculated energy spectra 
(general systematics), as well as charge radii, p 2 ( 0^ —> 0]*“), quadrupole moments for the lowest 2 + and 4 + states, and 
B(E 2\ 2 1 —> 0]*") (with k,l = 1, 2). C omparisons with existing data have also been carried out for 190 Po and 194 Po in 
Ref. fll.i ] and for 196 Po in Ref. |l!4l ]). Because their formulation is within a mean-field context, invoking the concept 
of nuclear shapes (spherical, oblate, prolate), it is not always possible to compare with the present IBM-CM approach 
in which the model space is restricted to Op-Oh and 2p-2h excitations across the Z=82 closed shell. A few general 
remarks are in place. Regarding the energy spectra the main conclusion of the calculations are the overestimation of 
the 0 2 energy (which is predicted to be oblate), though the slope of the variation with respect to A is well reproduced 
and, moreover, that the calculation is not able to reproduce the almost constant value of the 2f excitation energy, 
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FIG. 7: Experimental excitation energies and absolute B(E2) transition rates for selected states in 190 198 Po. 


giving rise to a steady increase of the energy as a fun ction of A. Another very interesting outcome of the calculation is 
the spectroscopic quadrupole moment (note that in 127l ] an equivalent magnitude is provided, f} 8 ) for the 2^ and 2^ 
states. The calculation predicts a prolate shape for the 2+ state while oblate for the 2^" state up to N=106 (A=190), 
interchanging at this point their character, becoming almost spherical (but slightly oblate) for N=114 (A=198) and 
onwards. In the case of A=190, the calculation exhibits the presence of a prolate and an oblate band structure, very 
much like the IBM-CM model results. The results for the charge radii give an overall correct trend. However, the 
upsloping trend in (r 2 ) (relative to the value for 210 Po), starting at N=112 (A=196) is not well described within the 
BMF description as compared with the IBM-CM. 


E. Wave function structure: from the unperturbed structure to configuration mixing 

We start our analysis with the structure of the configuration-mixed wave functions along the yrast levels, expressed 
as a function of the [N] and [N + 2] basis states, as given in eq. Q. In Fig. [TUI i and Fig. [TUb . we present the 
weight of the wave functions contained within the [A]-boson subspace, defined as the sum of the squared amplitudes 
w k (J,N) = | | 2 , for both the yrast states, k = 1, and the k = 2 states (the latter are indicated with a 

dashed line) for spins J n = 0 + ,2 + ,3 + ,4 + in panel (a) and J = 5 + ,6 + ,7 + ,8 + in panel (b). The results exhibit an 
interesting behavior: the yrast states show rather an intruder character for mass A=190 that quickly are changing into 
a regular one in the mass interval A=194-200. For the 0 + states, they are fully mixed for A=190 and subsequently, 
the 0+ state is changing into a pure regular one, while the character of the 0^" state is essentially purely intruder for 
A=200, later on changing into a pure regular character for A=204. The behavior for the 2 + and 4 + states is very 
similar. The first (2+ or 4 + ) state starts as a dominant intruder state for the lighter isotopes, steadily changing with 
an increase of its regular component, becoming a fully regular configuration for A=200 and onwards to the heavier 
masses. The second 2 + and 4 + state start with a regular component of ss 50% and ~ 30%, respectively, reaching a 
maximal value for A=194 (« 70%), then decreasing when moving to A=200, increasing up to reach a purely regular 
character by mass A=202. Indeed, the Oj, 2%, and 4 \ states exhibit a similar trend. The other states exhibit a more 
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FIG. 8: Theoretical excitation energies and absolute B(E2) transition rates for selected states in 190 198 Po. 


erratic behavior, mainly because of the many crossings in the unperturbed energy spectra of the regular and intruder 
states, although appearing at very similar energies. Finally, note that the results for A=208 should all correspond to 
pure regular states because the intruder states appear at a much higher energy at N=126. Because of the constraint 
put on the value of A w+2 (see the discussion in section ITV Cl) , the intruder states appear too low in the energy spectra. 
Moreover, because of the reduced number of bosons in 208 Po (Z=84, N=124), 2, only regular states up to spin J=4 
can be constructed within the IBM. Consequently, it is not possible to describe the higher-spin (J=8 and beyond) 
states showing up in the mass region 204 < A < 208 in a reasonable way within the IBM, in view of the specific 
high-spin broken-pair states appearing at ~ 1.5-2 MeV. 

Regarding the energy systematics of the intruder states, one expects a parabolic shape centered around N=104, as 
is the case of Hg and Pb. However, this is not the case for Po, neither for Pt isotopes. The reason for not observing 
the parabolic shape is the rather strong interaction between regular and intruder configurations and the crossing of 
the intruder and regular states in the ground state. Therefore, it is very enlightening to calculate the energy spectra as 
a function of A switching off the mixing part of the Hamiltonian, i.e. keeping w^’ N+2 = w^' N+2 = 0. These spectra 
are depicted in Fig. 1111 where we show the lowest two regular and the lowest two intruder states for different angular 
momenta. One observes a rather flat behavior of the energy for the regular states, still indicating a down-sloping 
tendency when moving towards the lighter Po isotopes. The energy of the intruder states is smoothly decreasing 
down to neutron mid-shell (N=106). In this region, the parabola become very flat (A=190-194). This results mainly 
from the smooth change of the Hamiltonian parameters when passing from isotope to isotope. A striking fact is the 
almost degeneracy of the unperturbed regular and intruder 0 + states for A=190-194, with the intruder configuration 
becoming the lowest one in the energy spectrum for A=190-192. The crossing of a regular and a intruder configuration 
with the same angular momentum has a strong influence on the regular component of the states resulting from the 
full IBM-CM calculation (see Fig. [TUI) inducing an interchange of character between the two states. In the situation 
of the J n = 2 + unperturbed energy spectrum, the closest approach happens at A=196 and in Fig.[TU]one notices that 
the 2+2 states interchange their character at this point. The same happens for the 4 + states at A=196-198. 

A most interesting decomposition of the wavefunction is obtained by first calculating the wavefunctions within the 
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FIG. 9: (Color online) Panel (a): Experimental and theoretical quadrupole moment for the 8)" state (note that the theoretical 
values for A=206 and 208 are not shown due to the reduced number of valence bosons). Panel (b): Gyromagnetic factor for 
the 8)*" state (experimental data and theoretical results). 


N subspace as 


n, jm)n 9 = E c '( JM ) 


(9) 


and likewise for the intruder (or N+2 subspace) as 

*(m, JM)t +2 = E C T(- N + 2)^((sd)f +2 ; JM) 


( 10 ) 


defining an “intermediate” basis [53], |54|. This generates a set of bands within the Op-Oh and 2p-2h subspaces, 
corresponding to the unpert urbed bands that are ex tracted in schematic two-level phenomenological model calculations 
(as discussed in references fl00 l[T37Ul3lfl4M52i ). and indeed correspond to the unperturbed energy levels depicted 
in Fig. EU 

The overlaps jv (2, JM \ k, JM) and N+ 2 (m, JM \ k, JM) can then be expressed as, 


n(1, JM | k, JM) = J2 a i(J', N)c\(J-, N), 


( 11 ) 


and 


N+2 (m,JM | k, JM) =Y J b k J (J\N+ 2)cJ{J]N+ 2), (12) 

3 

(see expressions © and (fTOll ). In Fig. [12] we show these overlaps, but squared, where we restrict ourselves to the first 
and second state (k = 1,2) with angular momentum J 7I '=0 + ,2 + ,3 + ,4 + ,5 + ,6 + ,7 + ,8 + , and give the overlaps with the 
lowest three bands within the regular (IV) and intruder (TV + 2) spaces (l = 1, 2, 3 and m = 1,2,3). Since these figures 
are given as a function of mass number, one obtains a graphical insight into the changing wave function content. In 
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FIG. 10: (Color online) Regular content of the two lowest-lying states for each J value (full lines with closed symbols correspond 
with the first state while dashed lines with open symbols correspond with the second state) resulting from the IBM-CM 
calculation, as presented in figure [3] 



A 

FIG. 11: (Color online) Energy spectra for the IBM-CM Hamiltonian presented in Table HT1 switching off the mixing term. The 
two lowest-lying regular states and the lowest-lying intruder state for each of the angular momenta are shown (full lines with 
closed symbols for the regular states while dashed lines with open symbols are used for the intruder ones). 


panel (a), which corresponds to the first state, one observes for the 0 + state a strong mixing between the first regular 
and first intruder unmixed states (also called the “intermediate basis”) for A=190-194, while for the heaviest isotopes 
the state is mainly of regular character. For the 2 + states the situation is very much the same though the mixing 
is much more reduced. For the higher spin states one observes a sharp border at A=198, in such a way that for A 
< 196 the states correspond to the first intruder configuration, while for A > 198 they correspond to the first regular 
configuration. Note that the intruder character observed for the heavier isotopes is somehow artificial because of the 
reduced number of bosons and also to the reduced excitation energy of the intruder configurations. In panel (b) the 
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Eh 1 reg. 
1 2 nd reg. 

B ~ rd 

3 reg. 
0 1 st int. 

2 nd int. 
B 3 rd int. 


(a) 



H 1 reg. 
1 2 nd reg. 

B 'i rd 

3 reg. 
0 1 st int. 

2 nd int. 
B 3 rd ini. 


(b) 


A 


FIG. 12: (Color online) Overlap of the wave functions of Eq. 0, with the wave functions describing the unperturbed basis 
Eq. (j9j) and Eq. GSJ. Panel (a): overlaps for first 0 + ,2 + ,3 + ,4 + ;5 + ,6 + ,7 + ,8 + state, panel (b): overlaps for the corresponding 
second state (see also text). 


situation becomes rather fuzzy especially for J > 2. For the 0 + state one observes a mixture between the first regular 
and intruder configurations up to A=196, then the composition changes into a pure first intruder configuration and 
finally ends up as the second regular configuration. Concerning the 2 + state, the lighter isotopes correspond to a 
mixture of the second intruder and the first and second regular configuration, the composition changing into the first 
intruder for A=198-200, finally ending as the second regular configuration. For the higher-spin states, the pattern 
is not so clear. In general, the states become rather pure and the changes in their composition is mainly due to the 
many crossings between the unperturbed regular and intruder configurations. 









































18 


V. STUDY OF OTHER OBSERVABLES: ALPHA-DECAY HINDRANCE FACTORS AND ISOTOPIC 

SHIFTS 

A. a-decay hindrance factors 

The a-decay process can be used as an important tool to probe nuclear structure, in particular the overlap factor 
between the initial and final wave functions has the character of an a- particle spectroscopic factor. Because calcula¬ 
tions of the absolute decay rate are very difficult, most often one studies decay branches, characterized by a hindrance 
factor described by the expression 


S 2 ip 

jjp _ 9 s _ 1 gs Jr ex 


s 2 


I P 

X ex ± gs 


(13) 


where 6 ? is the reduced a width, P ai the penetration probability through the combined Coulomb and centr ifuga l 
barrier |l06l | and /, the a-decay intensity (with i = gs, ex for the ground state and excited state, respectively) jl06| . 

The Pb-region has been studied intensively and hindrance factors have been used to obtain valuable information on 
shape coexistence, and mixing of various nuclear shapes. In particular, a two-level mixing model, suggested by Wauters 
et al. 97;], has been used to analyze the experimental hindrance factors and extract complementary information on 
mixing between different nuclear configurations. 


[N+2] 



[N+2] 



FIG. 13: Schematic view of a decay from the Rn nuclei into the Po nuclei, and from the Po nuclei into the Pb nuclei. 


The structure information on the composition of the 0+ ground state in the 19 °- 198 p 0 nuclei is derived from the 
specific a decay hindrance factors into both the 0+ and 0^" states into the corresponding daughter 186 - 194 Pb nuclei 
(see Fig. m- Use is made of the fact that for the Pb nuclei, it was shown, through study of the Pb nuclear charge radii 
@ and very weak E0 decay strength from the 0^~ -A 0[" transti on, that the Pb n uclei retain a spherical ground-state 
structure. The experimental results on these Po nuclei [95l - l97l . I99l - ll0ll Il06l Il07j have been analyzed in these papers 
resulting in hindrance factors of 2.8 ± 0.5, 2.5 ± 0.1, 1.1 ± 0.1, and 0.4 for A=198, 196, 194, and 192, respectively. 
The major result is the decreasing trend in these hindrance factors, even favoring decay to the excited 0^ excited 
state at 591 keV in 188 Pb starting from 192 Po. In those studies, and using a two -state mixing model [97] , an intruder 
character for the ground state 0[~ state in 192 Po turns out to be ss 63% [lOOl ] . or, even larger > 6 5% 101 1 . More 
recently, a-decay from 188 492p 0 h as b een observed into the final 0 + excited states in 184 > 188 pb ]lQ4 ], Their analyses 
results in hindrance factors into the first excited 0^ state of 0.08 ± 0.03 and 0.57 ± 0.12, respectively, which may 
well turn out as in indication of a very large overlap of the wavefunctions in initial and excited state in the final 
nucleus. In that latter paper, an analyses was carried out from inspecting the various energy minima in both the 
Po and Pb nuclei, obtaining qua litat ive information on overlaps of the corresponding wave functions, derived from 
the Nilsson-Strutinsky approach fl53 j|. An analysis of hindrance factors was also carried out using a similar method, 
making use of deformed mean-field wave functions at the energy minima of the total energy surfaces, pointing out 
that this approach, albeit different from an ap proach explicitely including mp-nh excitations across the closed Z=82 
shell closure, renders quite similar resuls Il54| . 

The wavefunctions, obtained from the IBM-CM (see also similar results on the study of the Pt fH] and Hg (64 ] 
nuclei on a-decay hindrance factors), give rise to the following theoretical values of 6%, 17%, 42%, 53%, corresponding 
with the weight factors of the intruder component (see Fig. m for the masses A=198, 196, 194, and 192, respectively. 
These numbers can be compared with the experimental relative hindrance factors as discussed before, with a rather 
good agreement. 
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On the other side, a-decay starting from the even-even A Iln nuclei (see Fig. into the A 4 Po nuclei, gives very 
strong fingerprints for observing excited 0 + states, in particular, experiments studying the decay of 202 Rn into 198 P o 
M and of 198 ’ 200 Rn[9i] into the 194 496 p 0 nuc i e i Th e analysis of these data, using a two-state mixing model [l48 |. 
resulted in detailed information on the more deformed intruder configuration in the 2f state. Its weight factor is 
changing from 7%, passing by 31% and 72%, towards the value of 88 %, for A=200, 198, 196, 194, respectively. It 
is interesting to compare these results with the present IBM-CM calculations resulting in a variation of the intruder 
weight factor of 5%, 15%, 47%, and 78% (see Fig. [TUI for the same mass values. This again indicates a rather good 
description of the changing wave function structure for the 2+ state in the Po even-even nuclei. 

Studies on a-decay from the odd-mass Po nuclei have been carried out f or the same mass region, with now an 
extra neutron coupled to the underlying Po even-even core [M [MEM [ml- These experiments show the presence 
of rather weakly particle-core coupled states, associated with a spherical configuration for the even-even underlying 
core, as well the appearance of a more strongly coupled structure, associated with coupling of the odd particle with 
the more deformed configuration of the even-even core. 




FIG. 14: (Color onli ne) P anel (a): Charge mean-square radii for the Po nuclei. Panel (b): Isotopic shift for the Po nuclei. The 
data are taken from 418]. 


A lot of effort has been invested into going beyond the more simplified two-state mixing model that has been 
extensively studied analysing the a-decay hindrance factors related to the decay into the excited 0 ~ st a tes in the 
daughter nuclei. In particular, more detailed calculations have been carried out by Delion et al. [155l 11561 ]. and 
more recently in [108J, emphasizing the need for a microscopic QRPA description that encompasses both nedutron 
and proton pairing vibrations as well as including proton 2p-2h excitations across the Z = 82 closed shell. More in 
particular, calculations have mainly concentrated in order to better understand the obs erved hindrance factors for 
a-decay leading into the neutron-deficient Po, Pb, Hg, and Pt nuclei. Karlgren et al. [l54| use the information on the 
energy minima in the (3 — 7 plane of the total energy surface as the essential input to calculate the hindrance factors. 
They also discuss a possible connection between the present deformed mean-field approach with a p icture in which 
specific particle-hole excitations are invoked within a shell-model approach. Moreover, Xu et al. 11571 ] formulates the 
a-decay hindrance factors within a density-dependent cluster model to describe the tunneling of a preformed cluster 
through the deformed Coulomb barrier. 
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A A 

FIG. 15: (Color online) Comparison of the value of /3 extracted from the theoretical quadrupole shape invariants (panel (a)) 
and the ones extracted from experimental and theoretical B(E2) values (panel (b)), corresponding to the 0 4 and O)!" states. 


B. Isotopic shifts 

Ex perimenta l information about ground-state charge radii is also available for both the even-even eh and odd- 
mass jlld Ill7| Po nuclei. Combined with similar data for the adjacent Po, Pb, and Pt nuclei, as well as for the 
odd-mass Bi, T1 and Au nuclei, the systematic variation of the charge radii supplies invaluable information on the 
ground-state wave function fiHMfini . We illustrate the overall behavior of ( r 2 ) A relative to the radius at mass A=210 
in Fig. [Till and the r elati ve changes defined as A (r 2 ) A = (r 2 )A+ 2 -{i~ 2 )a in Fig. ITTh. The experimental data are taken 
from Cocolios et al. [l 18| . 

To calculate the isotope shifts, we have used the standard IBM-CM expression for the nuclear radius 

r 2 = r 2 + PJjfrisrN + pNnppN + -P/v +2 (7iV+2-^ + (3 N+2 hd)pN+2 • (14) 

The four parameters appearing in this expression are fitted to the experimental data, corresponding to the radii of 
A=192-198 even-even Po isotopes |l!8j . The resulting values are 7 jv = —0.108 fm 2 , /3jv = 0.152 fm 2 , 7^+2 = —0.022 
fm 2 , and /3n+2 = 0.027 fm 2 . Note that in determining these parameters we took as reference point 204 Po instead of 
210 Po which is the reference for the experimental data. 

Panel (a) of Fig. [TT] shows the value of the radius relative to A=210. The first important fact is the strong 
deviation from the linear trend (marked with the dotted line) at A=198 an d downwards which is showing the onset 
of deformation confirmed in several recent experimental papers EH, EH. T° see this fad more clearly, we also 
present the value of the isotopic shift which enhances the appearance of the deformation. We mention the very good 
quantitative agreement between the IBM-CM and the experimental data, which confirms that the interplay between 
the [ N] and [N + 2] contributions in the 0 + ground state wave function along the whole chain of Po isotopes is well 
described by the model. Comparing with the data for the nearby Hg and Pt isotopes, one notices that the range of 
variation of the isotopic shift for Hg is only « 0.02 fm 2 , while for Pt this amounts to « 0.1 fm 2 and, finally, for Po 
becomes ~ 0.2 fm 2 . This fact is telling us how abrupt the onset of deformation in the Po is. Note the strong similarity 
of Fig. n~il> with Pt case [58j], where there is also an abrupt drop of the isotopic shift at A=184-186. 


VI. QUADRUPOLE INVARIANTS AND NUCLEAR DEFORMATION 

The IBM can provide us with both the energy spectra and the corresponding wave functions, as well as all derived 
observables (B(E2)’s, quadrupole moments, radii, ...), working within the laboratory frame, as well as the correspond- 
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FIG. 16: (Color online) Matrix coherent-state calculation for 190 ~ 194 Po, corresponding with the present IBM-CM Hamiltonian 
(table [II]). The energy spacing between adjacent contour lines equals 100 keV and the deepest energy minimum is set to zero, 
corresponding to the red color. 


ing mean-field energy surface, defining a nuclear shape over the /3 —7 intrinsic frame. This geometric interpretation of 
the IBM can be constructed using the intrinsic state formalism, as proposed in Refs. |l6ll - ll65 |. The extension of this 


regular and i ntruder c onfigurations was proposed by Frank et al, introducing 
1691. In Refs. 


formalism to describe simultane ously 

a matrix coherent-state method jl66l4l69| . In Refs. |59Ul70( a detailed description of the method and its application 
to Pt and Hg isotopes, respectively, can be found. In Figs. fl 6 h. 116b. and IT 6 b the IBM-CM mean-field energy surfaces 
are depicted for 190 Po, 192 Po, and 194 Po, respectively. In 190 Po the energy surface presents an oblate mini mum, that 
evolves into a very flat surface for 192 Po and in a spherical shape for 194 Po. Note that from Fig. 3 in Ref. jl27(, the 
192 Po and 194 Po nuclei exhibit an oblate minimum, however, in reaching 190 Po, the energy minimum switches over to 
become a prolate minimum, with a very similar deformation. 

Even though the shape of the nucleus is not an experimental observable, it is still possible to extract from the data 
direct information about various moments characterizing the nuclear shape corresponding with a given eigenstate. 
Using Coulomb excitation, it is possible to extract the most important diagonal and non-diagonal quadrupole and 
octupole matrix elements, including their relative signs and, in a model independent way, extract information about 
nuclear deformation as shown by Kumar, Cline and co-workers EzMzl- ,_, 

The essential point is the introduction of an “equivalent ellipsoid” view of a given nucleus [ml corresponding to 
a uniformly charged ellipsoid with the same charge, (r 2 ), and E2 moments as the nucleus characterized by a specific 
eigenstate fmllml. 

From the theoretical point of view the nuclear shape can be calculated using the quadrupole shape invariants. The 
quadrupole deformation corresponds to 


«M = V5(0,+ |[QxQ](°>|0+). (15) 

For the triaxial rigid rotor model Imi it is directly connected with the deformation parameter 

92 = q 2 , (16) 

where q denotes the nuclear intrinsic quadrupole moment. 

To calculate analytically the quadrupole shape invariants characterizing the nucleus in its ground-state and low-lying 
excited states, it is necessary to resort to a closure relation, 1 = Ylji m \JiM){JiM\, leading to the expression 

92, < = 5>+||Q||2+}<2+||Q||0+). (17) 

r 

In Table IIVI we present the theoretical value of q 2 corresponding to the Oj 1- and 0^~ states, for the whole chain of 
Po isotopes. In this case, only very few experimental values of the reduced E2 matrix elements (0^||Q||2+) have been 
measured |l!5| besides the strongest and dominating (O^HQI^ 4- ) matrix element. Because of the limited number of 
reduced E2 matrix elements, going into r=2,3,.. final states, it is not possible to extract reliable experimental values. 

For the same reason it is not possible to extract the triaxial shape variable <5, this is the reason for not including this 

variable in the discussion. The main conclusion of this table is that the deformation of both O^o states is very similar 
and is dropping when moving from mid-shell, at N=104, to the end of the shell at N=126 [115]. 
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TABLE IV: Quadrupole shape invariants for the 190 208 Po isotopes. 


Isotope 

State 

q 2 (e 2 b 2 ) 

Isotope 

State 

q 2 (e 2 b 2 ) 

190 

of 

6.2“ 

200 

of 

1.1“ 


0 2 + 

5.7° 


of 

0.9“ 

192 

Of 

oo 

p 

202 

Of 

1.1 “ 


Of 

4.5° 


Of 

0.9“ 

194 

Of 

3.6 

204 

Of 

0.7“ 


Of 

3.5 


Of 

0.7“ 

196 

Of 

1.7 

206 

Of 

0.6 “ 


Of 

2.8 


Of 

0.4“ 

198 

Of 

1.4 

208 

Of 

0.4“ 


Of 

1.2 


Of 

0.1“ 


“The effective charges have been taken from neighboring isotopes (see Table ns. 


_St arting fro m the quadrupole invariant (1151) . one can extract a value of the deformation f3 (see, e.g., references 

[MiMlIzi) 


3 Z e Tq A 2 / 3 ’ 


(18) 


where e is the proton charge and ro = 1.2 fm. Thus, we can extract values for f} corresponding to the ground-state 
Of and the first excited Of state. 

The resulting /3 values, extracted from (IT51) are shown in Fig. flST a'). and one notices an overall decrease, albeit with 
a plateau in the mass region 194 < A < 202, which is the region corresponding with an overall change in the energy 
scale for the first excited 2f state. 

We also provide a different measure of the deformation value and extract a value of /3 starting from a given B(E2) 
value through the expression: 


Q U7T y j —^ j — z) 

P ~ 3 Zerl A 2 / 3 (J020|J — 20}’ 


(19) 


where is a Clebsch-Gordan coefficient. 

In particular, we extract the value of (3, derived from the B(E2] 2f —> Of) value, but also for any other B(E 2; J + —> 
(J — 2)+) value along the yrast band. The extracted /3 values then give interesting information about a possible 
variation along the band. Note that we extract the value for Of from B(E 2; 2f —> Of) and for Of from B{E2\ 6f —y 
4f). The corresponding values are plotted in Fig. flSl bl and are fully consistent with the results, as shown in Fig. flSl al. 
They are slightly smaller because of using only a single but dominant B(E 2; 2 f —> Of) reduced transition probability. 

The results for the Po nuclei are contrasting with the nearby Hg and Pt isot opes. In a similar analysis as the 
one carried out here, but for the Hg isotopes (see Figs. 2 and 3 in reference [l80l ]'). the /3 values extracted from the 
IBM-CM results exhibit a strong difference in the interval 180 < A < 190. This is consistent with the /3 values 
extracted from the B(E2;2f —► O^ 1 ") (showing the less deformed oblate structure) and from the B(E 2; 6^ -A 4^) 
(associated with a more deformed prolate structure) and quite different from the overall smooth drop in the /3 values 
extracted in the Po isotopes. Comparing with the Pt nuclei, in which fewer data were available (see figure 3.27 in 1), 
down to mass A=188, j3 values extracted from the B(E2]2+ —> Of) and the B(E2;6f —> 4+) reduced E2 transition 
probabilities do not show the same large separation as in the case of the Hg isotopes. In the Pt case, a rather smooth 
transition (a bump) from the region with A > 188 moving towards the mid-shell region exists. This was shown to be 
a consequence of the rather strong mixing between regular and intruder configurations The behavior of 

within the context of the IBM-CM, indicated a similar bump as a function of mass number [591. 

Interesting information about the changing collective structure along the yrast band for the Po can be extracted 
studying the variation of the moment of inertia as one moves up the band s truct ure. A quantitative measure, that is 
often used, is the kinematic moment of inertia J^\ which is defined as MB 


J« = 


h 2 


dE 


d (t/ iyj A l)) 


h 2 (2J- l) 


2E 1 {J —>•(</ — 2)) ’ 


( 20 ) 
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FIG. 17: (Color online) Comparison of the experimental (a) and theoretical (b) kinematic moment of inertia for selected Po 
isotopes with 106 < N < 114. 




E y (MeV) E y (MeV) E y (MeV) E y (MeV) 


FIG. 18: (Color online) Comparison of the experimental moment of inertia in Po (a), Pb (b), Hg (c) and Pt (d), for selected 
isotones with 106 < N < 114. 


where E^(J —> ( J — 2)) is the energy diffence E(J) — E(J — 2). For an ideal rigid rotor, with a moment of inertia 
that is independent of the angular momentum J, the kinematic moment of inertia reduces to a constant. 

In constructing the experimental values, we use the yrast band energies for the corresponding J 7r values. 
However, approaching the high-spin levels with 8 + , 10 + ,12 + , it can happen that non-collective (mainly broken-pair 
configurations) do cross the collective structure, giving also rise to isomeric states as a fingerprint. Consequently, 
considering those states, the smooth variation of J*- 1 ) will change into a backbending pattern. Experimental data on 
lifetimes, y-ray intensities characterizing the decay within a band structure, as well the coincidence data with y-rays 
from the low-energy transitions in the experimental yrast ba nds, have been used to construct the resulting values. 
In particular, for the Po nuclei (see references [8i.l85H89l 90i . lllll . lll4f ). we cover the 106 < N < 114 (190 < A < 198), 
region, with the res ults prese nted in Fig. 1171 A similar analysis spanning a smaller set of Po nuclei, was carried out 
in references i mins, uni- One observes a very smooth increase along the yrast band, moving up to rather high 
spin values in the A=190 (up to spin 14 + ), 192 and 194 (up to spin 10 + ) isotopes, consistent with the energy spectra 
as shown in Fig. [l] It is only for the higher-spin states, which are of a non-collective broken-pair nature, that very 
small energy differences appear, not at all consistent with collective excitations and giving rise to strong back-bending 
(not shown here). From A=198, and onwards, a different energy pattern for J W is observed. For the lower part 
of the band, up to the 6^ level, the values of E 7 are roughly constant, indicating the typical energy differences of a 
vibrational structure. The fairly correct reproduction of this trend by the IBM-CM results indicates that the major 
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physics content is rather well described. 

It is also interesting to compare - making use of the experimental data on band structure and shape coexistence 
in the Pb region - the Po nuclei, with the neighboring Pb, Hg and Pt isotopes in order to appreciate similarities and 
differences, covering the same neutron interval 106 < N < 114, as shown in Fig. [TS] The data used to con struc t the 
jt 1 ' 1 band for the Po, Pb, Hg, and Pt isotopes have been taken from the references given in Section HT1 and jl82f . 

An obvious observation from Fig. [TH is the fact that at N=106, i.e. the nuclei 190 Po, 188 Pb, 186 Hg and 184 Pt, from 
J 7r = 6 + onwards, a very similar variation in the value of J W results. The starting point in these nuclei is ~ 30 h 2 
MeV -1 , followed by a similar slope up to spins 12 + and 14 + . In view of the well-established prolate character in the 
case of 188 Pb and 186 Hg, it is a convincing argument that the ground-state band in both 190 Po and 184 Pt, from spin 
4 + , 6 + , and onwards, behave like a prolate band too. One also notices that since the ground state 0 + and first excited 
2 + state in 188 Pb retain mainly a spherical character, the low E 7 part exhibits a different sloping start, very much like 
it is the case in 186 Hg. In the latter nucleus, there appears an off-set of the prolate band versus the ground-state less 
deformed band (oblate, anharmonic vibrational) causing almost constant values of J W to result for the 2 + ,4 + states. 
In 190 Po and 184 Pt, a very smooth behavior results along the yrast band, extending up to rather high-spin states and 
this for the whole region 106 < N < 114. One also observes a considerable drop in the value of in this interval, 
which is of the order of « 10 h 2 MeV -1 . In this respect, the Po and Pt nuclei exhibit a rather similar collective 
behavior, with the ’’intruder” states being not that obvious from the experimental data on excitation energies nor 
the corresponding B(E2) electromagnetic properties. Only radii and alpha-decay hindrance factors point towards the 
need to consider both configuration (regular [N] and intruder [N+2]) spaces. 


VII. CONCLUSIONS AND OUTLOOK 

The Po region constitutes one of the clearer examples where shape coexistence plays a key role in explainig the 
major observed features such as the evolving energy spectra, with mass number decreasing from the neutron closed 
shell at N=126, approaching the mid-shell region at N=104, the electromagnetic properties, in particular the B(E2) 
systematics, nuclear mean-square charge radii, a-decay hindrance factors, etc. These isotopes are situated in the 
region of the Pb, Hg and Pt isotopes where the proximity of the Z=82 proton closed shell, and its stabilizing effect 
to keep nuclei mainly into a spherical shape, in particular for the Pb nuclei, is very evident. These series of isotopes 
can be divided in two groups: on the one hand, the Pb and Hg nuclei, in which two (or even three for the Pb nuclei) 
types of configurations coexist and are experimentally well documented, characterized by energy spectra that exhibit 
the characteristic parabolic behavior of the intruder band structure, minimizing its excitation energy at the mid-shell, 
N=104 [52|, [59|. On the other hand, the Pt and Po isotopes for which one cannot disentangle easily the presence of 
two distinct band structures. Still, in both the Pt and Po isotopes, approaching the mid-shell point, the first excited 
0+ state behaves in an unexpected way, dropping seriously in excitation energy which can be used as an indirect hint 
for the presence of shape coexistence, however in a concealed way [58| ■ 

In the present paper, we have carried out an extensive study of the chain of isotopes 19 °- 208 p 0 using the Interacting 
Boson Model, including pair boson excitations across the Z=82 proton closed shell (called the intruder configuration 
space), and their interaction with the regular configuration space, called the IBM-CM. We have determined the 
Hamiltonian and the E2 operator describing this interacting system of bosons, through a least-squares fit to the 
known experimental data. This then results in the energy spectra, and, moreover, the calculation of many different 
observables such as the B(E2) values, nuclear mean-square charge radii, gyromagnetic factors and a-decay hindrance 
factors give the possibility to test the nuclear dynamics. In particular, the latter three observables are shown to 
serve as fingerprints to test the relative composition of the nuclear wave function into its ’’regular” and ’’intruder” 
components, and thus are strong indicators for the presence of shape coexisting structures. 

A very important issue concerns deriving information on the nuclear deformation properties. In particular, recent 
Coulomb excitation experiments on nuclei far from stability made it possible to extract a set of reduced E2 matrix 
elements (Of || T{E2) || 2j). These matrix elements can be used to construct the so-called ’’quadrupole shape 
invariants” for the 0+ and 0^" states. We have given particular attention to how /3-values can be derived from these 
invariants, or alternatively, from known experimental B(E2) values, and compared both approaches in the case of 
the Po nuclei. Thereby, a clear picture of the shape coexistence phenomenon in Po arises. There are two families of 
configurations, one slightly deformed (rather 7-unstable) or spherical that corresponds to the regular configurations, 
while the other more deformed and corresponding to the intruder configurations. Because of the lack of experimental 
information it is not possible to determine unambiguously the shape connected to the two families of configurations 
and, indeed, through the use of mean-field results (see section Hill) we assume an oblate shape for the intruder states, 
although the moment of inertia seems to suggest a prolate shape for the ground state of 190 Po. Both configurations 
have a rather large interplay and in the mid-shell region with the two 0 + unperturbed bandheads being degenerate 
in energy implying important mixing between the two configuration spaces. We have found that, in contrast with 
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the Hg nuclei, the difference in the quadrupole deformation extracted for the 0+ and 0^ states is quite small and the 
mass dependence very similar, independent of the method used to extract the quadrupole deformation. 

Concluding, it looks like the data available at present, based on our present study of the Po isotopes and a 
comparison with the Pt isotopes points towards rather similar structures. In both cases the ground-state at the 
mid-shell is composed mainly from the intruder configurations, although in the Po isotopes, only the second part of 
the neutron shell N=82-126 has been studied experimentally. For the Pt isotopes, almost the whole shell has been 
covered, but for the Po nuclei, because of the sparse character of the known data in and below N=106, A=190, a 
detailed comparison is very difficult. Moreover, the similar character of the Po and Pt isotopes is supported from a 
comparison of the kinematic moments of inertia for the Po and Pt isotopes: both exhibit a smooth variation as a 
function of neutron number in the interval 106 < N < 114. Comparing with the corresponding moments of inertia in 
the Pb and Hg isotopes, the nuclei 190 Po, 188 Pb, 186 Hg, and 184 Pt (starting at spin 6 + ) are very similar, resulting in 
the suggestion that from that spin, and onwards, the 190 Po ground band behaves like a prolate band. This is further 
corroborated by a comparison of the known experimental excitation energies for the prolate bands that have been 
observed in 186 Hg and in 188 Pb for states with spin from 6 + to 10 + . They match very well the corresponding energies 
in the yrast band structure in 190 Po. Possibilities to test this proposal might come partly from multiple Coulomb 
excitation with higher energy projectile ions (HIE-ISOLDE), and, possibly at a later stage, using one and two nucleon 
transfer reactions and studying E0 properties for the Po nuclei, which may lead to the direct observation of a low-lying 
excited 0 + state, as suggested by the present study. 
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